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Abstract 



We study the dynamics near the truncated p : ±q resonant Hamiltonian equilibrium for p, q 
coprime. The critical values of the momentum map of the Liouville integrable system are found. 
The three basic objects reduced period, rotation number, and non-trivial action for the leading 
order dynamics are computed in terms of complete hyperelliptic integrals. A relation between 
the three functions that can be interpreted as a decomposition of the rotation number into 
geometric and dynamic phase is found. Using this relation we show that the p : —q resonance 
has fractional monodromy. Finally we prove that near the origin of the 1 : —q resonance the 
twist vanishes. 
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1 Introduction 

The phase space of an integrable system is foliated into invariant tori almost everywhere. Exceptions 
occur when the integrals are not independent, i.e. when the energy-momentum map has critical 
values. Monodromy is a manifestation of the singular nature of the fibre over a critical value of 
the energy-momentum mapping. The term was introduced in Duistermaat's 1980 paper [7] as the 
simplest obstruction to the existence of global and smooth action-angle coordinates. Consider a 
closed, non-contractible loop T in the image of the energy-momentum mapping consisting entirely 
of regular values. Assume the loop is traversed from r(0) to r(l) where r(0) = r(l). For each 
regular value T{s) the period lattice Pr{s) gives the periods of the flows generated by energy and 
momentum, and thus encodes the transformation to action-angle variables. The period lattice Pr{s) 
at r(0) and at r(l) are related by a unimodular transformation. If this transformation is non-trivial 
the system is said to have monodromy. Cushman and Bates [6] use the ratios of periods, i.e. the 



1 



rotation number, to compute monodromy. Geometrically monodromy means that the torus bundle 
over the loop T is non-trivial. 

If the loop of regular values is contractible through regular values there is no monodromy. Thus 
there must be at least one critical value inside the loop T for monodromy to occur, and the nature 
of the critical fibre in principle determines the monodromy. For non-degenerate focus-focus points 
the critical fibre is a pinched torus with p pinches, and the monodromy matrix of a loop around 

the critical value of the focus-focus point is ^ |12| I20j . 

Recently Nekhoroshev et al. |15|[T6] found fractional monodromy in certain resonant oscillators. 
In fractional monodromy the change of basis of the period lattice is not an element of SL(2,7j) 
but of SL(2,Q) instead. This would be incompatible with the Liouville- Arnold theorem which 
gives the actions uniquely up to transformations in 5'L(2,Z). Thus for fractional monodromy the 
loop r necessarily has to cross critical values. The critical fibre has to be such that it is possible 
to continuously pass a sub-lattice of the period lattice through it [161 E] > a-lso see [10] . The first 
example was the 1 : —2 resonant Hamiltonian, with an appropriate compactification \15\ [TBI flT] - 
In [2] and [19] this has been extended to the compactified p : —q resonance. 

The present paper is based on the thesis [18]. It is a generalization from the 1 : —2 resonance as 
treated in [5] to the general p : —q resonance. We independently obtained results similar to those 
presented in jl9j . but our approach is complementary. There are two main differences. The first 
main difference is that we base our approach on the computation of the action variables. We derive 
the remarkable formula relating the action variables Ii and I2 by 

I^ = P±l^hT - hW (1.1) 
47r 

where A/i is the value of the non-linear part of the Hamiltonian, T is the period, and W is the 
corresponding rotation number. This formula can be interpreted as decomposing the rotation 
number into a dynamical phase proportional to T and a geometric phase proportional to the action 
I2 [13]. Then we show that W changes by l/(pg) upon traversing a loop around the equilibrium. 



and thus by (1.1) the action I2 changes by Ii / (pq) and the system has fractional monodromy. The 
second main difference is that we analyze the dynamics near the equilibrium point in detail. The 
period, rotation number, and action are dominated by certain singular contributions, and we show 
that sufficiently close to the equilibrium point these are the leading order contributions. This allows 
us to perform the computation without any compactification, thereby proving the conjecture made 
in [TB , that fractional monodromy is independent of the details of the compactification. Finally 
we obtain the new result that the isoenergetic non-degeneracy condition of the KAM theorem is 
violated near the 1 : —q resonance, i.e. the system has vanishing twist near the equilibrium point. 

The paper is structured as follows. First we derive the resonant Hamiltonian normal form 
near an elliptic-elliptic equilibrium point in p : —q resonance which is then studied in the reduced 
phase space after applying singular reduction. After discussing the structure of the reduced phase 
space we derive the set of critical values of the energy-momentum mapping. The period of the 
reduced flow, the rotation number, and the non-trivial action of the system are computed in the 
next section. We then prove fractional monodromy for the p : —q resonance. Finally we prove that 
the twist vanishes in 1 : — gr resonant systems close to the equilibrium point. 



2 Resonant Hamiltonain Equilibria 

Near a non-resonant elliptic equilibrium point a Hamiltonian can be transformed to Birkhoff normal 
form up to arbitrary high order. The truncated Birkhoff normal form depends on the actions only. 
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and to lowest order these actions are those of the two harmonic oscillators of the diagonalized 
quadratic part of the Hamiltonian. When the frequencies are p : zizq resonant, the normal form 
contains additional so called resonant terms that depend on resonant linear combinations of the 
angles. Even though the non-linear non-resonant terms are generically present, we will assume that 
all the non-linear non-resonant terms up to the order p+q are vanishing Without this assumption the 
dynamics near the equilibrium would be dominated by the non-resonant terms when p + q > 5. For 
the low order resonant cases no such assumption is necessary, and therefore the results presented 
in [5] for the 1 : —2 resonance are relevant for generic Hamiltonian systems. By contrast, the 
integrable systems studied in this paper have higher codimension (increasing with p + q) in the 
class of Hamiltonian systems. 

Denote a point in phase space T*M^ with coordinates {pi,qi,P2,Q2) and let the symplectic 
structure be w = dqi A dpi + dq2 A dp2- Assume the origin {pi, qi,P2, 92) = (0, 0, 0, 0) G T*M? is an 
elliptic equilibrium point of the system whose eigenvalues +iuji, +iuj2 have ratio p/q, where p and 
q coprime positive integers. Then the quadratic part of the Hamiltonian H near the equilibrium 
point can be brought into the form 



by scaling uji to p by a linear change of time. For a = +1 the quadratic Hamiltonian H2 is definite, 
while for o" = — 1 it is indefinite. The system is said to be in p : +q resonance, and mostly we are 
interested in the case a = —1. 

The classical treatment of resonant equilibria, see e.g. [2], is as follows. Denote {Ai,(f)i), i = 
1,2, the canonical polar coordinates corresponding to {pi,qi)- The resonant Birkhoff normal form 
depends on Ai, A2 and the resonant combination (the so called secular term) —aqcpi + p4>2- The 
normal form Hamiltonian truncated at order p + q is 



where the phase (p can be set to zero by a shift of the angles. In order to reduce the number of 
variables a linear symplectic transformation to {Ji,ipi) is performed. The resonant combination is 
introduced as a new angle ^2, while the cyclic angle ipi is conjugate to Ji = pAi + aqA2 = H2. The 
complete transformation is given by J = M~^A, tp = M(j), where M contains two arbitrary integers 
a and b restricted by det M = bp — aaq = 1 . The new Hamiltonian depends on Ji, J2, ip2 only, 
and is therefore integrable. The lowest order term in 'i/'2 has coefficients proportional to y^^A^A^. 
Setting the non-linear non-resonant terms /i^j ec[ua.l to zero the Hamiltonian becomes 

H = Ji + /iA/(6Ji - aqJ2)''{-aJi + pJ2)P cos'ilj2 ■ 
Reverting back to the original Euclidean coordinates the resonant term reads 



and the Hamiltonian simply is H = H2 + fiAH. The functions H2 and AH are in involution and 
independent almost everywhere, and thus the system is Liouville integrable. This is the integrable 
system we are going to analyze in this paper. Note that unlike previous work |16 |llll[T9] . we do not 
add higher order terms to AH to compactify the system. Moreover, we will study the integrable 
system defined by the Hamiltonian H, instead of the Hamiltonian AH. For the discussion of 
monodromy the two are equivalent, since monodromy is a feature of the Liouville foliation, that 
does not depend on the dynamics on the individual tori. Also, they share the same singularly 





H = pAi + aqA2 + fJ^ijA{Al + ^lyj A\Al cos(-ct#i + p(t)2 + 'p) 



AH = ^[{pi+iqiy{p2-aiq2Y] 
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reduced system. However, we are also interested in the physically relevant dynamics of the Birkhoff 
normal form in the context of KAM theory, and therefore we analyse the Hamiltonian H and not 
the Hamiltonian AH. In particular, when considering the rotation number of the full system the 
difference is crucial. Nevertheless, we will study the momentum map F = {AH, H2) (which is not 
the energy- momentum map for the Hamiltonian H), and the Hamiltonian H = H2 + fiAH is a 
(linear) function of the momenta. 



3 Reduction 

In this section we review the steps necessary to reduce to a single degree of freedom using the 
symmetry i?2- We state only the results and refer the interested reader to the standard literature, 
for example Cushman [6], Abraham and Marsden [1] and Broer et. al. |4j, for the p : —q resonance 
in particular also see fTD]. These steps retrace the derivation reviewed in the previous section in a 
more geometric way. 

The flow of the resonant quadratic part H2 is the S'^-group action 

: X — >C^ 

{t,zi,Z2) ^ {ziexp{pit) , Z2exp{aqit)) , zi, Z2 e C 

where Zi = -^{pi+ iqi), i = 1, 2. This action is non-degenerate except on the axis zi = and 

Z2 = on which points have isotropy subgroup Zg C and Zp C 5^, respectively (if q,p > 1). 
The invariants of this group action are 

VTi = ZiZi, 7r2 = 22^2, ^3 = ^^(-^l-^f)' = ^(z^z^) (3.3) 

for cj = — 1. For a = +1 instead vrs = ^(2:^22)) ^4 = ^(-^1^2) invariants. 

Because of the non-trivial isotropy we employ singular reduction in order to reduce the system 
to one degree of freedom. The reduced phase space 

Ph, = H,\h2)/S' 

is given by the relation 

ttI + ttI = nfTT^, vri > 0, 7r2 > 0, (3.4) 

hence it is a semi-algebraic variety in the invariants. The ambient Poisson space is endowed with 
the Poisson structure 

{7ri,7r2}=0, {vri,7r3} = g7r4 

{vri,7r4} = -gvra, {vr2, vrs} = -cjp7r4 

{tt2,7T4} = pairs, {vrs, 7r4} = l'^i''^'^2~^ (crp^vri - g^7r2) 



with Casimir 7i2 = pui + q(JTT2 and symplectic leaves the reduced phase spaces P} 



h2- 



Fixing the Casimir 7^2 = h2 we choose to eliminate 112 from (3.4) and thus get an equation for 
the reduced phase space in the form 

,2 , _2 _ q fh2-p7riY 
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Figure 1: The section 7r4 = of the reduced phase space P/12 ^or the 2 : ±3 resonance. Upper left: 
a = +1 and /12 = 1.1. Upper right: a = —1, /12 = 0.5. Lower left: a = —1, /12 = —0.5. Lower 
right: o" = — 1, /12 = 0. For /i2 7^ points vri = are cusp singularities, and points vri = /i2/2 are 
conical singularities, both due to the non-trivial isotropy of 
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The interval of valid vri is determined by the requirements tti > and 7r2 > 0, thus 
o- = +1 : TTi G [0, h2/p] 

a = —1, /i2 > : TTi G [/i2/p, oo) (3.6) 
(J = -1, < : TTi e [0, oo). 

Fig. [T] shows sections 714 = of the rotationally symmetric reduced phase space P^^ for the 
2 : ±3 resonance for all four relevant cases 0" = +!, 0" = — 1 and /12 > 0, /12 < and /i2 = 0. For 
a = +1 the reduced phase space is compact with a cusp-singularity {q = 3) at the origin vri = 
and a conical singularity (p = 2) at vri = /12/2 (where 7r2 = 0). These singularities in the reduced 
phase space are a result of the non-trivial isotropy of the group action <I>^2 at zi = vri = and 

Z2 = TT2 = 0. 

For cr = — 1 the reduced phase space is non-compact. The singular points are of the same type 
as in the compact case, but they exist separately for positive or negative /i2. For /i2 > there is a 
conical singularity {p = 2) at 7r2 = and for /12 < there is a cusp singularity [q = 3) at tti = 0. 
For /i2 there is a singularity of order p + g = 5 at vri = (in the compact case a = 1 the reduced 
phase space is merely a point for /i2 = 0) . 



It follows from equation (3.5) that in the general p : ibg resonance there is a singularity of order 
q at TTi = and of order p at 112 = 0, assuming /12 7^ 0. If /12 = the order of the singularity is 
p + q. Note that in the above statement a singularity of order 1 means no singularity. 

Expressed in the invariant the integral AH simply becomes tt^, so that the reduced Hamiltonian 

is 

^(vri, 7r2, VTs) = pm + aqTT2 + /ivrs = /12 + /Uvrs . (3.7) 

As mentioned earlier, the truncated resonant Birkhoff normal form of a generic resonant Hamilto- 
nian system would also contain terms vr^^vrg, with i + j < {p + 9)/2, but in order to maximize the 
effect of the resonant term these are assumed to be zero. 



4 The Energy— Momentum Mapping 

Let 

F : {pi,qi,P2,q2) ^ {H2{pi,qi,P2,q2),^H{pi,qi,p2,q2)) 

be the momentum map of the integrable system, and denoted its value by (/i2, A/i). The elliptic 
equilibrium at the origin is a critical point of the momentum map F since both integrals are of 
order > 2 in {pi,qi), so that lank DF{0) = 0. Furthermore, this point is degenerate in the sense of 
the momentum map [3]. For this we need to show that the Hessians of H2 and AH are not linear 
independent at a; = 0. This follows immediately from the fact that AH is of order p + q > 3 in 
{pi,qi) so that its Hessian vanishes identically at the origin. 

The critical values (/i2,A/i) of the momentum mapping for the p : zizq resonance (i.e. the 
bifurcation diagram) are shown in fig. [2| and are described in the following lemma. 

Lemma 1. The entire line Ah = is critical for p > 1, q > 2 and a = —1. For p = 1, the critical 
values are given by Ah = and /i2 < and for a = +1 they are given by Ah = and /12 > 0. 
Moreover, there are two additional branches emanating from the origin of the bifurcation diagram 
for a = 
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Figure 2: The critical values (/i2, Ah) of the energy-momentum mapping are shown as red branches. 
The top left picture corresponds to the 1 : —q resonance, the top right one to the p : —q resonance 
with p > 1. The degenerate equilibrium point is marked by a disk at the origin for a = —1. The 
lower left picture corresponds to the 2 : 3 resonance with a = +1. 



Proof. A point in reduced phase space is critical if 1) it is a singular point of the reduced phase 
space or 2) if it is a point of tangency of the surface ATC = Ah (which is simply a horizontal plane) 
and the reduced phase space. The first condition together with (3.6) gives that tti = vra = 7r4 = 
is critical for (j/i2 > and that tti = h2/p,7r3 = 774 = is critical for /12 > 0. The corresponding 
critical values of the momentum map are (/i2)0). The second condition requires the gradients of 
the two integrals to be parallel, thus 714 = 0, vrs 7^ 0, and the derivative of the right hand side of 
(3.5), vr^7r2, with respect to vri vanishes. This gives p^vri = aq'^iT2. Since vri,7r2 are non-negative 
there is no tangency for a = —1 (appart from the singular point with /i2 = 0). For a = 1 there are 
two additional families of critical values emerging from the origin in a cusp of order p + q. □ 



5 Dynamics near the Degenerate Equilibrium 

In this section we derive equations for the period T of the reduced flow, the rotation number W 
of the full system, and the second action I2. A prominent feature of the reduced period T is its 
algebraic (rather than logarithmic) divergence when approaching the degenerate equilibrium point. 
This is especially easy to see once we introduce weighted polar coordinates {p, 9) in the bifurcation 
diagram. In these coordinates, T and W separate into p- and ^-dependent contributions which 
considerably simplify the computations. 
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5.1 The Reduced Period 



An equation for T is derived by separation of variables from 

■fri = {TTi,7i} = iJ.{TTi,ir3} = fiqir^. (5i 



Using the equation for the reduced phase space (3.5) together with Ah = /^tts, tt^ can be written 



as a polynomial in vri. It follows that the reduced period is defined on the hyperelliptic curve 

r = {(7ri,u;) G I = Q(vri)} where Q{z) = fi^'^-^z'' {a{h2 - pz)f - {qAhf 
Separating the variables in equation ( |5.8[ ) and integrating yields 

T{h2,Ah)= <f^. (5.9) 
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Our main focus is the p : —q resonance with non-compact fibration. Thus the integral along a 
closed loop as it stands makes no sense. We define the reduced period by dividing the dynamics into 
two parts: Dynamics close to the equilibrium point, and dynamics far away from it. If the system is 
compactified by an appropriate higher order term (as in |16| [TT] ) we may assume that the dynamics 
far away from the equilibrium will eventually return to the neighbourhood of the equilibrium. The 
time spent on this return loop is a smooth function of initial conditions if we assume that there 
are no additional critical points in it. Specifically we consider a Poincare section at some small 
but finite value of vri intersecting stable and unstable invariant manifolds, similarly to the analysis 
done for symplectic invariants in |17| [9] and the 1 : — 2 resonance in [5]. The contribution of the 
near-dynamics (from the section of the stable manifold to the section of the unstable manifold) 
is divergent when approaching the equilibrium point, while the contribution of the far-dynamics 
(from the section of the unstable manifold to the section of the stable manifold) remains a smooth 
and bounded function. 

For convenience we modify the truncated period thus defined one more time. Notice that 
the integral of dTTi/w from any finite positive tti to oo is finite, and smoothly depends on the 
parameters. Thus changing the truncated period integral to an integral over the non-compact 
domain only changes it by a smooth function, and the same argument applies. As a result we can 
treat the closed loop integral of the non-compact system as our leading order period. In particular 
it correctly describes the leading order divergent terms when approaching the equilibrium point. 

An alternative point of view that combines the last two steps (first restriction to the near- 
dynamics, then the extension to oo) is to consider the integral in a compactification of the complex 
plane into a projective space. In this space the integral for T is compact, and this also explains 
why it is bounded in the first place. This approach was first used in [5] . 

In [IH] two approaches are followed. The first one generalises the treatment of the 1 : —2 
resonance [HI [TO] to the p : —q resonance. There is a privileged compactification which prevents the 
rotation number from diverging when approaching the critical values. Using this compactification 
the period lattice, i.e. reduced period and rotation number, is computed. Then the authors use the 
Newton polygon to find the leading order terms of these integrals in the limit approaching the origin 
of the bifurcation diagram. The analogue of our polynomial Q appears in their work as the Newton 
polygon approximation. What our approach shows is that the resulting hyperelliptic integrals 
have direct dynamical meaning as explained above. In the second approach in [19j the problem is 
complexified and monodromy is found as Gauss-Manin monodromy of a loop with complex (/i. Ah) 
that avoids critical values. As the authors point out the singularity at the origin is not of Morse 
type, and this is related to the fact the singularity is degenerate in the sense of the momentum 
map. As a result the period diverges algebraically (instead of as usually logarithmically). 
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Lemma 2. The reduced period T diverges algebraically with exponent | A/i| p+i upon approaching 
the degenerate equilihrium point on a curve with non-vanishing derivative at the origin. For p = 1, 

the period diverges like \h2\ ^ on the line Ah = when /i2 ^ 0. 

Proof. Consider the polynomial Q as a polynomial in tti, /12, and Ah. It is weighted homogeneous, 
where /i2 and tti must have the same weight, so that their weight is 2 while that of Ah is p + q. 
Therefore we introduce weighted polar coordinates {p, 9) in the image of the momentum map by 

A/i = /+«sin0 (5.10a) 
h2 = p^cose (5.10b) 

Together with vri = p^x it follows that 

Tih2, Ah) = p-(P+^)+^ / — =: p-(P+«)+2A(e) (5.11) 

where 

= p^q^-Px'i {a{cose - px))P - {qsinOf 

is independent of p. Thus T factors into a radial and an angular contribution A(9). The transfor- 
mation for (/i2) Ah) to (p, 9) is at the origin but C°° everywhere else. 

The period T diverges with p'^~P~^. When approaching the origin on a line with non- vanishing 
slope (in the original variables (/12, A/i) the contribution in Ah ~ fP^'^ is of leading order. Thus 
T ~ I A/ip/*^^'^'^)^"'^. When p = 1 the line segment Ah = with /12 > is not critical. Approaching 
the origin along this line there is no contribution from Ah, and thus only /i2 ~ is relevant, and 
the result follows. □ 

5.2 The Rotation Number 

Recall from the introduction the symplectic coordinate system with {iJ;i,Ji} = 1. The reduced 
period gives the period of 1^2, while the rotation number gives the advance of ^/'i/(27r) during that 
period. The ODE for ipi is 

i^i = {9i,H}. (5.12) 

Integration of this ODE for time T gives the rotation number. Note that the angle "01 is not be 
globally defined, but all we need is the derivative of ipi, see the comments in [5J. 

The interpretation of the rotation number is similar to the interpretation of the period: The 
true (compactified) rotation number will differ from ii by a smooth function. The leading order 
singular part of the rotation number is contained in R. 

Proposition 3. The rotation number R of the p : zizq resonance is given by 

.,..A.).^/(l..f(^-^))^^, 

Proof. The angle tpi satisfies the following Poisson brackets: 
= b 

{V'i,vr2} = -a 

vra / qb ap^ 



{V'ljTTs} 



TTl 7r2 
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These brackets follow from expressing tTj in terms of the canonical variables ( Jj, ipi), noting vri = Ai, 
-7r2 = A2, TTs = ^^''^^2''^ cos ^^2, and then computing the canonical brackets, using 

diTi dlT2 



Thus, the differential equation for tpi (5.12) becomes 
TTs / qb ap^ 



2 V TTi TT2 



Changing the integration variable from t to tti using (5.8) gives Tpi as an Abelian integral. Com- 
parison of the limiting behaviour of the rotation number and insisting on the relations R = uoxjuj^ 
and dH/dIi = cui gives the overall sign —a in R. □ 

We refer to the three integrals R is composed of as 

R{h2,Ah) = -a (^T{h2,Ah) + W{h2, Ah) 
y zvr 

W{h2,Ah) = Wiih2, Ah) + W2ih2, Ah) 

where 

W,{h2,Ah)= ^Ah^'^^' 



47r / -Kiw 
W2{h2,Ah) = -'^Ah^'^''' 



4vr / 'K2W 

Expressing these functions in the weighted polar coordinates (p, 0) gives 

B^{e) = ^sme I — (5.13) 
47r J xw 

B,(,),-gsin0/ , (5.14) 

47r J [px — cos Ujw 

Note that Bi and B2 do not depend on p, which is the main virtue of the weighted polar coordinates. 
Nevertheless, the original functions Wi and W2 are not continuous at (/i2, Ah) = (0, 0). The reason 
is that they take different values when approaching the origin along different lines 9 = const. 

Figure [3] shows a plot of Bi{9) on the left and B2{0) on the right on a loop T{6) around the 
origin of the bifurcation diagram fig. |2]top right with constant p for the 2 : —3 resonance. 

The functions Bi and B2 are periodic, but discontinuous. The discontinuity occurs where the 
critical points are crossed. The functions can be made continuous by an appropriate shift, but then 
they will not be periodic any more. They cannot be made periodic and continuous at the same 
time. It turns out that this behaviour is the reason fractional monodromy occurs in the p : —q 
resonance, see section [6j From the derivation of Bi and B2 it is clear that they are not really 
independent functions, but can be obtained from each other by exchanging p and q, and shifting 9 
by IT. 

Note that if we would have considered AH as our Hamiltonian instead of H the rotation 
number would not have the diverging contribution from T; it would he W = Wi + W2 instead of 
R = T/2'K + W . For the vanishing twist described later on this contribution is crucial. 
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Figure 3: Plot of Bi{9) (left) and B2{0) (right) for the 2 : -3 resonance, a = -1, 6 = 2. The 
corresponding bifurcation diagram is shown in fig. [2] top right. Bi is discontinuous on the branch 
of critical values with Ah = and h2 < 0, i.e. 9 = itvr, with total jump equal to — b/q. B2 is 
discontinuous on the branch of critical values with A/i = and /12 > 0, i.e. = 0, with total jump 
equal to — a/p, see Lemma [5] 



5.3 The Non-Trivial Action 

The flow of H2 is periodic with period 27r, and therefore Ii = H2 is one action of the system. The 
other action I2 is a non-trivial function of A/i (or h) and Ii. From the local canonical coordinates 
mentioned in the introduction we know that J2 and tjj2 are conjugate variables. Translating this to 
the reduced system gives 

Proposition 4. The non-trivial action I2 of the p : ibg resonance is 

hiluAh) = -a (^^^AhT{h,Ah)- hW{h,Ah)^ . (5.15) 

Proof. Using the reduced Poisson bracket we find 
{cos~^ —^^=,aTTi + hTT2} = -a . 



TT. 
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The action I2 is therefore given by 

l2(/i2, Ah) = (h arccos .^^ (advri + 6d7r2). 
2vr J V^i7r2 

Integration by parts gives 



l2 = — i n^^ Z/aq al.r\ ^,^+^&P + ^-2 



27r 



Now TTi and t^2 are related by the Casimir 1\ = h2 = pvri + aq'K2- Thus d'K2 = a^dni and with 
Ah = yUTTs this equation becomes 

—a /Ah, , , f diTi Ahh2ap f dni Ahh2bq f dni 

h = -7r{v + q) {bp - aaq) f + f f 

27r \ 2 J w 2 J 7r2W 2 J ttiw 

Using det M = bp — aaq = 1, /i2 = /i, and recalling that W = W1 + W2 gives the result. Notice that 
the final answer is independent of the choice of integers a and b in M, as long as det M = 1. □ 
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This expression for I2 reduces to the one found for the 1 : —2 resonance [5], where p + q = 3, 
a = and 6=1. In the next section on fractional monodromy we shall prove that the terms Wi 
and W2 are discontinuous (see fig. [3| at = ibvr and 9 = ±0 respectively. They can be made 
non-periodic and continuous, thus causing the action I2 to be globally multivalued. 
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Even though T diverges algebraically like Ahp+i ^ (Lemma 2), the action I2 goes to zero like 

p"^ ~ Ahp+1 when approaching the equilibrium point. The action I2 does have the interpretation of 
a phase space volume. Even though the system is non-compact the action does not diverge when 
approaching the equilibrium point. The geometric reason is that I2 measures the volume relative 
to the (unbounded) separatrix, which is finite. The boundary terms from the partial integration 
cancel, see [18] for the details. 

Another interpretation of this formula is obtained by solving it for the rotation number. This 
gives a decomposition of the rotation number into a dynamical phase proportional to T and a 
geometric phase proportional to the action I2, compare e.g. [13] . 



6 Fractional Monodromy 

In this section we establish the fact that the p : —q resonance has fractional monodromy. We 
explicitly calculate the monodromy matrix M that gives the transformation of the actions Ii and 
I2 after one full anticlockwise cycle around a loop T enclosing the degenerate equilibrium point at the 
origin of the bifurcation diagram fig. [2] top row. We then describe the singular fibres corresponding 
to the critical values the loop T crosses. 

Lemma 5. The function Bi{9) satisfies 

lim BiiO) = ± — . (6.16) 

The function B2{9) satisfies 

limB2{e) = T^. (6.17) 
e^o± 2p 



This is the main technical result, but instead of duplicating the proof here, we refer to [19] . 
The main addition to their work at this point is the interpretation of these integrals. In [19] they 
appeared as the leading order Newton-polygon approximation of the compactified rotation number 
integrals. In our approach they appear directly as the integrals of the rotation number of the 
non-compact system interpreted as explained before. 

Using this technical result we are now giving another proof of fractional monodromy in the 
p : —q resonance which is based on the explicit expression of the action obtained earlier. Recall 
that Ii is the globally smooth action H2 of the system as introduced in the introduction. 

Theorem 6. The p : —q resonance has fractional monodromy near the degenerate elliptic equilib- 
rium point. The actions change according to 




after one full anticlockwise cycle on a loop T around the degenerate equilibrium point in the 
bifurcation diagram. 
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Proof. Assume the loop T with fixed p = po'is traversed in the mathematical positive sense. When 
crossing the line A/i = 0, /12 < 0, "^^^^ ^ = and ^^^2-^ ^ = Hence, the effective jump 

of Wi becomes — |. A similar argument holds for W2 at A/i = 0, /12 > 0. Although the loop is 
traversed such that 9 crosses the line A/i = 0, /i2 > from below, W2 has the opposite sign of Wi. 
Thus, the effective jump of W2 is — ^. Using the form of the action variable from proposition |4j it 
follows that the action I2 changes like 

T ( 0. h 
l2^h-h\- + - 



upon completing a full cycle. Due to 

a b aq + hp 1 
p q pq pq'' 

the monodromy is (independently of the integers a and b) given by 

I2^l2--h- 

pq 

□ 

The singular fibre of the critical values with /i2 < depends on the value of q. Consider the 
Poincare section P2 = and p2 < 0. The critical level has Ah = vrs = 0. Thus we need to 
study the level set determined by the three equations pj-zip + crg|z2p = /i2 < 0, ^^(z^Zg) = 0, 
and P2 = near the critical point zi = 0. From the first and last equation we obtain ^2 = 
(— /i2 + p\zi\'^)/{—aq). For small \zi\ therefore q2, and hence Z2, is approximately constant. The 
remaining equation thus becomes ^{zf) = 0. Writing zi = rexp'np gives rf cos qip = 0, so that the 
level sets of the intersection of the critical fibre with the Poincare section are given by the 2q rays 
with qip = 7r(n + 1/2), n = 0, 1, . . . , 2g — 1. Together they form q straight lines which are the stable 
and unstable manifolds in alternation, all passing through the origin. 

If the level set is compactified, e.g. by considering ^{zf) + |zi|''"'"^, it becomes a flower with q 
petals. The action of the flow of H2 on the level set is as follows. Start with a point in the section 
P2 = 0, hence Z2 = q2- It returns to the section when the imaginary part of its image vanishes 
again with positive derivative, hence for Q{q2exp{aqit)) = q2sm{aqt) = 0, with aq q2COs{aqt) < 0. 
The smallest positive t that solves this is t = 27: /q. The action of the flow after the return time 
therefore is 

(2vr/g, 2:1,^2) = (ziexp(2i7rp/g),Z2) . 

This is simply a rotation by 2TTp/q in the zi plane. As a result the petals of the flower are mapped 
into each other, and p is the number of petals that the rotation advances by. Since p and q are 
coprime, all petals are visited before the orbit returns back to the initial one. The action is the same 
for either sign of a, the difference is that for a = —1 instead of petals we have sectors delineated 
by stable and unstable manifolds. 

Considering the level set as a whole, and not only in the Poincare section, gives a curled torus 
whose transversal cross section is the petal, and whose tubes rotate by 2'Kp/q when they complete 
one longitudinal cycle. The fact that we discussed the flow of H2 instead of the flow of H (or /S.H) 
is immaterial for the topology of the level set since the flows commute. 

For /i2 > in a similar way a flower with p petals appears in the section pi = 0. 
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All critical values are degenerate when q > 2. This is a crucial difference between the 1 : — 2 
resonance and the higher order resonances. In particular the terms mj that were set to zero in 
the normal form could completely change the bifurcation diagram. How to envisage the singularity 
at the equilibrium point with /12 = A/i = is unclear. Both curled tori limit to this unstable 
degenerate equilibrium point, but the above argument breaks down since |/i2| cannot be assumed 
large compared to \zi\ any more. 



7 Vanishing Twist in the 1 : —q Resonance 

For the p : —q resonance with p > 1 the bifurcation diagram is divided into two halves by the 
critical line Ah = 0. We suspect that in this case the twist does not vanish for regular values, 
but we have not been able to find a proof of this. However, when p = 1 the line of critical values 
stops at the origin, see fig. |2j This is therefore a typical case where one can expect vanishing twist 
to occur, see ^ for a general topological proof. In this particular case using the weighted polar 
coordinates allows for a simple analytical proof of vanishing twist. 

Theorem 7. For the 1 : —q resonance, q>2, the twist vanishes near the degenerate equihbrium 
point on the curve Ah = 0, /i2 > 0. 

Proof. By definition the twist vanishes if the rotation number R has a critical point on the energy 
surface, i.e. when R does not change from torus to torus. On the energy surface this is written as 
dR/dIi\ 

H=const- In the image of the energy momentum map the condition is satisfied when the 
lines of constant rotation number and the lines of constant energy are tangent to each other. This 
implies that Vi? and 'VH are parallel. The gradients can be computed in any coordinate system, 
and we choose the weighted polar coordinates. Thus the condition for vanishing twist is 

dR dH dRdH _^ 
dp do do dp 

Now H = H2 + p.AH where H2 is order 2 m. p and AH is order p + q in p. For small p we can 
therefore neglect AH, and find dH/dp ~ 2/9 cos 0, dH/dO ~ — p^sin^. Recalling the factorisation 
T = p-^P+i^+'^A{0) from ^1a\ we find 



d fi" dT^ 

-^T27r^ = 1^ = {-{p + q) + 2)p~^P+^)+^A{0), 

-cr27r^ = % + B'AO)+ B'M ^ p-^^^'^^^^ ^ {0) . 
The leading order contribution comes from T alone. Altogether this gives 

up to lowest order in p. The solution ^ = follows from symmetry: A{0) is even in thus A'{0) is 
odd in 0, and thus both terms are odd in 0. For p = 1 the line = is a line of regular values, and 
therefore we have shown that the twist vanishes along ^ = asymptotically near the origin. □ 
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